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Abstract. In the present paper we prove that every k-dimensional Cartesian 
product of Kingman convolutions can be embedded into a k-dimensional sym- 
metric convolution (k=l, 2, . . . ) and obtain an analogue of the Cramer-Levy 
theorem for multi-dimensional Rayleigh distributions. 
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1. Introduction, Notations and Preliminaries 

In probability theory and statistics, the Rayleigh distribution is a continuous 
probability distribution which is widely used to model events that occur in differ- 
ent fields such as medicine, social and natural sciences. A multivariate Rayleigh 
distribution is the probability distribution of a vector of norms of random Gauss- 
ian vectors. The purpose of this paper, is to introduce and study the fractional 
indexes multivariate Rayleigh distributions via the Cartesian product of Kingman 
convolutions and, in particular, to prove an analogue of the Levy-Cramer theorem 
for multivariate Rayleigh distributions. 

Let V := denote the set of all probability measures (p.m.'s) on the positive 

half-line K + . Put, for each continuous bounded function f on M + , 

(1) [°° f(x)fi *i j v{dx) - [{S " " 



o v^r( s + i) 

OO />OG nl 

/ / f((x 2 + 2uxy + y 2 ) 1 / 2 )(l-u 2 y~ 1 / 2 n(dx)v(dy)du, 

JQ J-l 

where \i and v £ V and 5 — 2(s + 1) > 1 (cf. Kingman [7i and Urbanik [IT]). The 
convolution algebra (V, *i,s) is the most important example of Urbanik convolution 
algebras (cf Urbanik [17] ). In language of the Urbanik convolution algebras, the 
characteristic measure, say a s , of the Kingman convolution has the Rayleigh density 

(2) da 8 (y) = 2(8 + * )S+ V s+1 exp (-(s + l)y 2 )dy 

F(s+1) 

with the characteristic exponent x = 2 and the kernel A s 

(3) A s (x) =r( S + l)J s (x)/(l/2a;) s , 
where J s (x) denotes the Bessel function of the first kind, 

{-l) k {x/2f+ 2k 



(4) J„(x) := S^o- 



k\T{v + k + 1) 
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It is known (cf. Kingman [7], Theorem 1), that the kernel A s itself is an ordinary 
characteristic function (ch.f.) of a symmetric p.m., say F s , defined on the interval 
[-1,1]. Thus, if 9 S denotes a random variable (r.v.) with distribution F s then for 
each teR+, 



(5) A s (t) =Eexp(it9 s ) = J cos (tx)dF s (x). 

Suppose that X is a nonnegative r.v. with distribution and X is independent 

of 6 S . The radial characteristic junction (rad.ch.f.) of /i, denoted by £i(t), is defined 

by 

/• DC. 

(6) jl(t) = Eexp(itXd s ) = / A s (tx)iJ,(dx), 

Jo 

for every t G K + . The characteristic measure of the Kingman convolution *i t g, 
denoted by o~ s , has the Maxwell density function 

,„s daJx) 2(s + l) s+1 o,,-. , , „. o, ,„ , 

(7) —fa = r( s + i) x ex Pi-( s + 1 ) x }> (0<x<^). 

and the rad.ch.f. 

(8) a s (t) =exp{~t 2 /4(s + l)}. 

Let P := V(M. + ) denote the class of symmetric p.m.'s on R + . Putting, for every 
GeP, 



F S (G)= / F cs G(dc), 
Jo 



we get a continuous homeomorphism from the Kingman convolution algebra ("P, *i,s) 
onto the ordinary convolution algebra (V, *) such that 

(9) F s {G 1 * 1 ,sG 2 } = (F S G 1 )*(F S G 2 ) G 1 ,G 2 eV 

(10) F s a s = N(0,2s + 1) 

which shows that every Kingman convolution can be embedded into the ordinary 
convolution *. 



2. Cartesian product of Kingman convolutions 

Denote by R +fc ,fc = 1,2,... the k-dimensional nonnegative cone of M. k and 
V(M. +k ) the class of all p.m.'s on K +fc equipped with the weak convergence. In the 
sequel, we will denote the multidimensional vectors and random vectors (r.vec.'s) 
and their distributions by bold face letters. 

For each point z of any set A let S z denote the Dirac measure (the unit mass) at 
the point z. In particular, if x = (xi, x 2 , ■ ■ ■ ,Xk) € R fc+ , then 

(11) <5 X = S Xl x 8 X2 x . . . x 5 Xk , (k times), 

where the sign " x " denotes the Cartesian product of measures. We put, for 
x = (xi, ■■■ ,x k ) and y = (2/1,2/2, • • • ,Vk) € K +fc , 



(12) S x Q s , k 5 y = {5 Xl o s S yi } x {S X2 o s 5 V2 } x ■■ ■ x{S Xk o s S Vk }, (k times), 
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here and somewhere below for the sake of simplicity we denote the Kingman con- 
volution operation *i,s,S = 2(s + 1) > 1 simply by o s ,s > 5. Since convex com- 
binations of p.m.'s of the form (fTTj) are dense in T'(M. +k ) the relation (fT2")l can be 
extended to arbitrary p.m. 's Gi and G 2 € T(R +k ). Namely, we put 

(13) G 1 Q s ,fe G 2 = JJs x Qs, k 5 y G 1 (dx)G 2 (dy) 



which means that for each continuous bounded function cf> defined on R +fe 
(14) I #a)Gi Q s ,k G 2 (dz) = ^ { I cf ) (z)S x O s . k S y (dz)}G 1 (d x )G 2 (d y ). 

R+fc R+ fc R+ fc 

In the sequel, the binary operation Os.fc wm be called the k-times Cartesian prod- 
uct of Kingman convolutions and the pair (T'(R +k ), Os,fc) will be called i/ie fc- 
dimensional Kingman convolution algebra. It is easy to show that the binary op- 
eration Q s .k is continuous in the weak topology which together with ([T]) and (fTB")) 
implies the following theorem. 

Theorem 1. The pair (T'(M. +k ), Os,fe) * s a commutative topological semigroup with 
So as the unit element. Moreover, the operation Qs,fc is distributive w.r.t. convex 
combinations of p.m.'s inV(M. +k ). 



For every G € V(R +k ) the k-dimensional rad.ch.f. G(t),t = t 2 , • ' ■ *fc) e 
° +fc , is defined by 



A- 



(15) G(t) = | nA s fex,)G(dx), 

R+fc 3 =1 

where x = (21, x 2 , ■ • • x&) € ffiL +fe . Let S = {9 s ,i,6 s ,2, • • ■ , where 6* SJ - are inde- 
pendent r.v.'s with the same distribution F s . Next, assume that X = {Xi, X 2 , X^} 
is a k-dimensional r.vec. with distribution G and X is independent of & s . We put 

(16) [0 S , X] - {B a ,iX l ,e s , 2 X2, s , k X k }. 
Then, the following formula is equivalent to (p~5|) (cf. [13] ) 

(17) G(t)=Ee i<t + & " x ->, teR +fe . 

The Reader is referred to Corollary 2.1, Theorems 2.3 & 2.4 [T3] for the principal 
properties of the above rad.ch.f. Given s > —1/2 define a map F s ^ k : "P(lR +fc ) — > 
P(E fc ) by 

(18) F s , fc (G) = /" (T Cl F a ) x (T C2 F S ) x . . . x (T Cfc F s )G(dc), 

R+ fe 

here and in the sequel, for a distribution G of a r.vec. X and a real number c we 
denote by T C G the distribution of cX. Let us denote by V s , k (R +k ) the sub-class of 
"P(R fc ) consisted of all p.m.'s defined by the right-hand side of (fT5|) . By virtue of 
(fl"5)) - (fl"8")) it is easy to prove the following theorem. 

Theorem 2. The set f> s , k (R +k ) is closed w.r.t. the weak convergence and the 
ordinary convolution * and the following equation holds 

(19) G(t) = 7-(i^, fe (G))(t) teR +k 
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where .F(K) denotes the ordinary characteristic function (Fourier transform) of a 
p.m. K. Therefore, for any Gi and G2 G M. +k 

(20) F.,n(Gi) * F fl , fc (G 2 ) = F a , fc (Gi Oa.k G 2 ) 

and i/ie map i^fc commutes with convex combinations of distributions and scale 
changes T Cl c > 0. Moreover, 

(21) F s>fe (£ 5 , fe ) = iV(0,2(s + l)I) 

where £ Sj fc denotes the k- dimensional Rayleigh distribution and N(0, 2(s + 1)1) is 
ffte symmetric normal distribution on R fe u/ii/i variance operator R — 2(s + 1)1,1 
oemg i/ie identity operator. Consequently, Every Kingman convolution algebra 
(7 , (M +fe ), Os,fc) * s embedded in the ordinary convolution algebra (l : ' sk (M. +k ), -k) 
and the map F Si k stands for a homeomorphism. 

Proof. First we prove the equation (|19[) by taking the Fourier transform of the 
right-hand side of (fig)) . We have, for t g R fc , 

J-(F s , fc (G))(t) = /n^cosfexjO^^G)^ 

(22) = / / cos(^)(T Cj .F s (dx)G(dc) 



J Y[A s (t jCj )G(dc) 



3 = 1 

= G(t) 

which implies that the set set 'P Si / c (R +fe ) is closed w.r.t. the weak convergence and 



the ordinary convolution * and, moreover the equations (|20j) and ([211) hold. □ 

Definition 1. A p.m. F g V(M. +k ) is called Q)s,k — infinitely divisible (Q s ,k—ID), 
if for every m—1, 2, . . . there exists F m £ P(K +fe ) such that 

(23) F = F m Os,fc F m Os,fc ■ • • Os,fc F m (m times). 

Definition 2. F is called stable, if for any positive numbers a and b there exists a 
positive number c such that 

(24) T a F Cj s , k T b F = T C F 

By virtue of Theorem [2] it follows that the following theorem holds. 

Theorem 3. A p.m. G is Q s ,k —ID, resp. stable if and only if H s .k(G) is ID, 
resp. stable, in the usual sense. 

The following lemma will be used in the representation of Os.fc — ID, k > 2. 

Lemma 1. (i) For every t > 

(25) hm = lim 1^!! = t. 

(ii) For any x = (xq, x\, ■ ■ ■ , Xk) and t = (to, t\, ■ ■ ■ , £&) £ K fe+1 , k = 1, 2, ... 

(26) Zzm^ 1 ~ n *^ Afl( ' ra:r) = l -^X{Arg{*))tl, 
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where p — ||x||, Arg(x.) — t^ui an d X r (Arg(~x.)),r = 0, 1, k are given by 
(27) \ r (Arg(x)) 



cos</> r = 0, 

sin (f> sin 4>\ ■ ■ ■ sin </> r _i cos 4> r 1 < r < k — 2, 

sin0sin0i... sin(/)fc_2 cos?/' f = k — 1, 

v sin sin (/>!... sin </)fc_ 2 sin ^ r = k, 



where < tp, <p, (f> r < 7r/2, r = 1, 2, k — 2 are angles ofn appearing its polar form. 

The following theorem gives a representation of rad.ch.f.'s of Os,fe— ID distribu- 
tions (see [13| ), Theorem 2.6 for the proof). 

Theorem 4. ^4 p.m. /i G /D(Os,fc) */ a^rf onfo/ i/ t/iere exist a a-finite measure 
M (a Levy's measure) on R +fc with the property that M(0) = 0, M is finite outside 
every neighborhood of and 

(28) / ^_M(dx) < oo 

J R+k 1 + ||x|| 2 

and for each t = (t\, tk) G R +fe 



(29) 



logA(t) = / {1 -JJUtjXj)} 1 *}^ 2 M(dx), 

JR+ k - =1 ll x ll 



where, at the origin 0, the integrand on the right-hand side of V29]) is assumed to 
be 

(30) ^U X J t2 J = Km|M,_ {l - 11 A s(t^)}^S~ 



j = l " " 

for nonnegative Xj,j — 1, 2, k given by equations in Lemma{^ In particular, 
if M = 0, then p becomes a Rayleighian distribution with the rad.ch.f (see definition 

w 



(31) -\ogm = \Y. x ^l teR +fc , 

3=1 

for some nonnegative Xj,j = 1, k. Moreover, the representation i29\) is unique. 

An immediate consequence of the above theorem is the following: 

Corollary 1. Each distribution p € /-D(Os.fc) * s uniquely determined by the 
pair [M,A], where M is a Levy's measure on R +fe such that M(0) = 0, M is 
finite outsite every neighbourhood of and the condition i28\) is satisfied and 
A := {Ai, A2, • • • Afe} G R +fc is a vector of nonnegative numbers appearing in 1131)) . 
Consequently, one can write p = [M,A]. 

In particular, if M is zero measure then p = [A] becomes a Rayleighian p.m. on 
R +fe as defined as follows: 

Definition 3. A k-dimensional distribution, say E s ,fc, is called a Rayleigh distri- 
bution, if 

(32) S s ,j; = cr s x a s x • • • x a s (k times). 



6 



THU VAN NGUYEN 



Further, a distribution F G "P(R +fc ) is called a Rayleighian distribution if there exist 
nonnegative numbers X r , r = 1, 2 • • • k such that 

(33) F = {T Xl a s } x {T X2 a s } x . . . x {T Xk a s }. 

It is evident that every Rayleigh distribution is a Rayleighian distribution. More- 
over, every Rayleighian distribution is Os,fc~ ID. By virtue of ([7]) and (|32[) it follows 
that the k-dimensional Rayleigh density is given by 

(34) g(x)=Ilj =1 ^ s + i) x? +1 exp{-(s + l)M 2 }, 

where x = (x\, xi, • ■ • > £fc) € R +fe and the corresponding rad.ch.f. is given by 

(35) £ Sik (t)=Exp(-\t\ 2 /4(s + l)), tel+ fc . 
Finally, the rad.ch.f. of a Rayleighian distribution F on M. +k is given by 

(36) F(t) = Exp(-±Y, X jfy 

3=1 

where Xj, j = 1,2, . . . , k are some nonnegative numbers. 

3. An analogue of the Levy-Cramer Theorem in multi-dimensional 
Kingman convolution algebras 

We say that a distribution F on M. k has dimension m, 1 < m < k, if m is the 
dimension of the smallest hyper-plane which contains the support of F. The follow- 
ing theorem can be regarded as a version of the Levy-Cramer Theorem for multi- 
dimensional Kingman convolution. The case k=l was proved by Urbanik (|18j). 

Theorem 5. Suppose that G 4 e V(R +k ),i = 1,2 and 

(37) = Gi Os.fc G2. 

Then, Gj,i = 1,2 are both Rayleighian distributions fufilling the condition that 
there exist nonnegative numbers \i >ri i = 1, 2 and r = 1, 2, . . . , k such that for each 
i=l, 2 the number of non-zero coefficients \ >r s among A^i, A,,2, • ■ • , \,k are equal 
to the dimension of Gi, respectively. Moreover, 

(38) AL + 7 L = 1, r = l,2,...,fc 
and 

(39) Gi = T Kl a s x T\ i 2 o~ s x . . . x T Xi k <j s 

Proof. Suppose that the equation (|37p holds. Using the map F St k we have 

F s , fe (S S;fc ) =F a , fc (Gi)*F aifc (G 2 ) 

which, by virtue of (|2T|) . implies that 

iV(0, 2(a + 1)1) = F s , fe (G!) * F s , fe (G 2 ). 

By the well-known Levy-Cramer Theorem on R k (cf. Linnik and Ostrovskii [9]), 
that they are both symmetric Gaussian distributions on R fc . Consequently, they 
must be of the form (l3"9l and the coefficients Aj, r ' 's satisfy the above stated condi- 
tions. □ 
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